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ABSTRACT 

An  optically  pumped  Cesium  vapor  magnetometer  was  used  to 
measure  the  temporal  fluctuations  of  the  total  geomagnetic 
field  at  Monterey,  California.   Power  spectra  were  obtained 
from  these  recordings  for  the  frequency  range  0.1-10   Hz. 
Measurements  were  made  during  the  month  of  April  19  78  with  the 
Kp  index  varying  from  0+  -  8+.   The  power  spectra  obtained 
displayed  a  characteristic  -20  dB/decade  slope  from  0.1  -  2.0 
Hz  and  were  essentially  flat  between  2.0  -  10.  Hz.   Measurement 
times  were:   Local  night  (0000-0200),  local  morning  (0800-1000) 
and  local  afternoon  (1600-1800).   Local  morning  exhibited  the 
highest  power  density  levels,  exceeding  the  lowest  power  den- 
sity levels  (recorded  during  local  night)  by  lOdB  at  .1  Hz  and 
3dB  at  10  Hz.   A  comparison  of  a  magnetically  disturbed  day 
(Fredricksburg  a  index  of  7)  with  a  magnetically  quiet  day 
(Fredricksburg  a  index  of  3)  at  the  same  local  time  (1730-1930) 
showed  that  at  the  low  frequency  end  of  the  spectrum  (.1  -  .6 
Hz)  there  was  lOdB  more  power  on  the  disturbed  day. 
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I.   INTRODUCTION 

It  is  often  convenient  to  consider  the  geomagnetic  field 
as  made  up  of  two  major  components:   (1)  The  main  field  which 
can  be  considered  to  be  approximately  dipole  and  generated  in 
the  interior  of  the  earth,  and  whose  observable  temporal  fluc- 
tuations are  on  the  order  of  years  or  longer.   For  purposes  of 
this  study,  we  shall  consider  the  main  geomagnetic  field  to  be 
constant  in  time.   (2)  A  time-varying  component  superimposed 
on  the  main  geomagnetic  field  whose  origin  is  to  be  found  in 
the  ionosphere  and  the  magnetosphere.   This  time-varying  com- 
ponent of  the  geomagnetic  field  is  of  small  amplitude  compared 
to  the  main  field  (1  part  in  10 *  or  less)  and  involves  frequency 
components  up  to  the  kHz  range. 

The  study  of  this  time-varying  component  of  the  earth's 
magnetic  field  is  of  interest  in  gaining  a  better  understanding 
of  basic  geophysical  processes  as  well  as  in  improving  the  per- 
formance of  existing  or  projected  devices  which  are  employed  by 
the  Navy.   A  number  of  operational  Navy  systems  such  as 
M(agnetic)  A(nomaly)  D(etection) ,  magnetic  mines,  torpedoes  and 
others,  as  well  as  some  projected  systems  (such  as  ELF  Communi- 
cation Systems),  depend  on  the  detection  of  small  changes  in 
the  existing  magnetic  field  in  the  ocean.   Thus,  from  both  a 
scientific  and  a  Naval  applications  point  of  view,  detailed 
experimental  data  on  the  geomagnetic  field  fluctuations  seem 
highly  desirable. 


As  a  first  step  in  an  experimental  program  to  characterize 
and  understand  more  fully  the  magnetic  field  fluctuations  in 
and  near  the  ocean,  we  have  measured  the  power  spectrum  of 
the  fluctuations  in  Monterey,  California  using  an  optically 
pumped,  total  field  magnetometer.  -The  measurements  were 
carried  out  at  a  relatively  quiet  field  site  near  La  Mesa 
village.   A  total  of  5  2  hours  of  geomagnetic  field  recordings 
were  made,  at  various  times  of  day  and  under  varying  conditions 
of  geomagnetic  activity.   Frequency  analysis  was  performed  on 
the  data  yielding  the  power  spectrum  of  the  fluctuations  in 
the  frequency  range  of  0.1  to  10.0  Hz.   The  general  shape  of 
the  power  spectra  agrees  with  earlier  measurements  but  several 
novel  features  were  also  observed. 
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II.   BACKGROUND 

A.   SOURCES  OF  GEOMAGNETIC  FLUCTUATIONS 

Figure  1  shows  an  approximate  power  spectrum  of  the  fluc- 
tuations measured  on  the  earth's  surface  and  having  periods  of 
a  day  or  shorter.   Also  indicated  on  the  graph  are  the  approxi- 
mate amplitudes  of  the  field  fluctuations  showing  that  even  the 
largest  of  these  fluctuations  are  less  than  1  percent  of  the 
magnitude  of  the  Main  Field  (B  -  50,000nT).   Nevertheless,  they 
are  of  considerable  interest  both  from  the  point  of  view  of 
geophysics  and  from  the  point  of  view  of  Navy  applications, 
and  a  world-wide  network  of  magnetic  observatories  maintains 
continuous  records  of  measurements  of  these  varying  field 
components.   World  Data  Center  A  (operated  by  NOAA  in  Boulder, 
Colorado)  issues  several  indices  as  well  as  forecasts  of  geo- 
magnetic activity  on  a  daily  basis. 

The  major  diurnal  variations  have  a  period  of  24  hours  and 
are  caused  by  current  systems  generated  in  the  upper  atmosphere 
as  a  result  of  solar  heating.   At  mid-latitudes  this  Sq  effect 
tends  to  be  2  5  -  50  nT,  but  in  the  equatorial  regions,  the 
effect  of  the  equatorial  electro jet  can  cause  a  typical  vari- 
ation of  about  one  hundred  nT  even  during  magnetically  quiet 
conditions.   There  also  exists  a  much  smaller  diurnal  variation 
with  a  period  of  12  hours,  caused  by  lunar  atmospheric  tides. 
The  magnitude  of  this  variation  is  less  than  0 . 1  of  the  daily 
solar  variation. 
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Figure   1.       Power  spectrum  of  geomagnetic  disturbances  observed  on  the 
Earth's  surface. 

[Cladis,    (1977)] 
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Spherical  harmonic  analysis  of  these  so-called  Sq  and  L 
variations  from  a  world  network  of  observatories  allows  one  to 
determine  the  equivalent  electric  current  systems  responsible 
for  the  observed  magnetic  variations.   About  two-thirds  of  the 
variations  are  found  to  be  due  to  currents  flowing  external  to 
the  earth,  and  the  remaining  third  can  be  attributed  to  currents 
induced  in  the  earth's  surface  layers  by  the  variable  external 
field.   The  external  currents  are  found  to  flow  mostly  at  an 
altitude  of  about  100  km  and  are  produced  by  convective  move- 
ments of  the  conducting  upper  atmosphere  across  the  earth's 
magnetic  field  lines.   By  analogy  with  the  electric  dynamo, 
this  mechanism  is  often  referred  to  as  the  atmospheric  dynamo. 

Magnetic  storms,  which  are  severe,  world-wide  magnetic 
disturbances,  have  their  main  frequency  components  in  the  same 
region  of  the  spectrum  as  diurnal  variations.   However,  the 
rate  of  occurrence  of  magnetic  storms  seems  to  be  primarily 
influenced  by  the  solar  sunspot  cycle. 

Magnetograms  obtained  at  low  and  mid-latitude  stations 
indicate  that  many  storms  have  a  characteristic  pattern  of 
development  as  the  storm  progresses.   The  start  of  a  typical 
storm,  signalled  by  an  abrupt  world-wide  increase  in  the  hori- 
zontal component  of  the  magnetic  field,  is  known  as  a  sudden 
commencement.   This  increase  is  typically  20  -  30  nT  and  has 
a  rise  time  of  a  few  minutes.   The  next  2-8  hours,  during 
which  the  field  remains  at  values  higher  than  pre-storm,  is 
known  as  the  initial  phase.   This  is  followed  by  a  12  -  24  hour 
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period  during  which  the  field  decreases  to  values  typically 
100  nT  below  its  undisturbed  magnitude,  and  this  is  known  as 
the  main  phase.   The  final  stage  of  the  storm  is  characterized 
by  a  gradual  return  of  the  field  to  its  pre-storm  value  and 
lasts  one  to  three  days. 

Individual  storm  records  will,  of  course,  show  deviations 
from  this  "typical  storm",  particularly  at  higher  latitudes, 
where  rapid,  large-scale  fluctuations  are  common  during  storms. 
Even  at  mid-latitudes  large  amplitude  fluctuations  with  periods 
of  about  one-half  hour  are  occasionally  observed  during  the 
main  phase  of  a  storm. 

Sudden  changes  in  the  dynamic  pressure  of  the  solar  wind 
impinging  on  the  outer  boundary  of  the  geomagnetic  field  are 
believed  to  be  responsible  for  the  initiation  of  magnetic 
storms  and  the  subsequent  development  represents  the  response 
of  the  magnetosphere  and  the  charged  particles  within  it  to 
this  initial  impulse. 

The  frequency  region  of  greatest  interest  to  us  lies  below 
3  Hz,  commonly  known  as  ULF.   Geomagnetic  field  fluctuations 
in  this  frequency  region  are  usually  referred  to  as  micro- 
pulsations,  with  periods  varying  from  0.2  -  600  seconds,  and 
amplitudes  from  a  fraction  of  an  nT  at  the  high  frequencies  to 
several  tens  of  nT  at  the  low  frequency  end  of  the  interval. 
Micropulsations  are  divided  into  two  general  types :   Continuous 
or  pc  and  irregular  or  pi.   The  pc  micropulsations  display 
amplitude  variations  which  are  quasi-sinusoidal  and  last 
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several  hours,  whereas  pi  micropulsations  have  considerable 
irregularities  both  in  frequency  and  amplitude.   As  indicated 
in  Figure  1,  the  spectrum  of  micropulsations  has  been  sub- 
divided into  roughly  5  segments  ranging  from  pc  1  at  the  high 
frequency  end  to  pc  5  at  the  lowest  frequencies .   The  occur- 
rence of  pc ' s  has  a  diurnal  variation  with  a  maximum  near  local 
noon.   It  is  believed  that  pc's  represent  small  hydromagnetic 
disturbances  existing  in  the  daytime  magnetosphere,  but  the 
possibility  that  they  are  of  ionospheric  origin  cannot  be  ruled 
out  at  this  time.   The  pi's  are  trains  of  micropulsations  con- 
sisting of  series  of  damped  oscillations  lasting  from  a  few 
minutes  to  about  one  hour  and  having  periods  of  roughly  40  -  100 
seconds.   Their  amplitude  is  largest  in  the  auroral  zone  and 
decreases  with  decreasing  latitude,  and  they  appear  to  be  pre- 
ferentially generated  in  the  night-time  magnetosphere.   At  high 
latitudes,  there  is  evidence  that  pi's  are  hydromagnetic  waves 
which  have  been  greatly  amplified  by  plasma  instabilities 
created  in  the  magnetosphere  by  high-energy  electrons. 

In  the  ELF  frequency  range  (3  -  3000  Hz),  a  number  of 
emissions  contribute  to  the  magnetic  fields,  the  most  pronounced 
ones  being  the  earth-ionosphere  cavity  resonances  (Schumann 
resonances).   These  signals  result  from  disturbances  that  are 
resonantly  excited  by  lightning  transients  in  the  concentric 
spherical  cavity  formed  by  the  earth's  surface  and  the  lower 
region  of  the  ionosphere.   The  power  spectra  of  the  ELF  signals 
often  show  maxima  near  7.8,  14.1,  20.3,  26. 4  and  3  2.5  Hz 
[Schumann  and  Konig,  19  54]. 
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B.   REVIEW  OF  EARLIER  WORK 

Power  spectrum  analysis  of  the  geomagnetic  field  has  been 
a  popular  tool  to  better  understand  geomagnetic  field  activity 
and  to  gain  insight  into  possible  source  mechanisms  for  the 
fluctuations.   For  many  naval  applications  one  must  also  under- 
stand these  geomagnetic  phenomena,  so  that  in  attempting  to 
detect  man-made  magnetic  anomalies,  one  might  try  to  eliminate 
this  background  signal.   This  allows  for  concentration  on  a 
smaller  frequency  range,  which  is  of  importance  to  the  military 
observer.   The  next  few  pages  are  devoted  to  a  short  review  of 
work  done  in  power  spectrum  analysis  of  the  geomagnetic  field. 

In  1963,  Santirocco  and  Parker  reported  [Santirocco  and 
Parker   1961]  their  work  on  the  polarization  and  power  spectra 
of  pc  micropulsations .   Recordings  of  the  N-S ,  E-W,  and  verti- 
cal components  of  the  geomagnetic  field  in  the  frequency  band 
.00  5  -  5  Hz  were  made  at  a  shore  location  in  Bermuda.   The  field 
sensors  utilized  were  ferromagnetic-cored  loop  antennas ,  polari- 
zed so  that  dB/dt  vectors  pointing  North,  East  and  down  were 


Geomagnetic  micropulsations  are  field  variations  with 
periods  ranging  from  .2  seconds  to  approximately  10  minutes. 
A  terminology  recommended  by  the  International  Association  of 
Geomagnetism  and  Aeronomy  is  pc  for  pulsations  continuous  and 
pi  for  pulsations  irregular.   The  period  ranges  are: 


pc 

1 

.2 

-   5 

sec 

pi  1 

1 

-  40  sec 

pc 

2 

5 

-  10 

sec 

pi  2 

40 

-150  sec 

pc 

3 

10 

-  45 

sec 

pc 

4 

45 

-150 

sec 

pc 

5 

150 

-600 

sec 

For  further  details  see  Jacobs  [1970 ]  ,  p.  16-20. 
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defined  as  positive.   The  sensitivity  of  the  loops  was 
70[uv-sec/nT]  (where  1  nT  =  1  Gamma  or  1CT5  Gauss).   With  350 
hours  of  continuous  taping  of  the  field  components  reported, 
power  spectrum  analysis  of  the  tapes  was  carried  out.   The 
power  spectra  of  the  pc  pulsations  appeared  as  clumps  of 
variable  fine  structure  components  superimposed  on  a  residual 
background  which  falls  monotonically  with  a  slope  of  -6dB/ 
octave.   Santirocco  and  Parker  also  observed  that  upon  com- 
parison with  the  results  of  others  that  this  background  slope 
is  constant  from  .000  2  Hz  to  nearly  1  Hz. 

Davidson  reported  in  1964  results  on  the  average  diurnal 
characteristics  of  geomagnetic  power  spectra  in  the  period 
range  4.5  -  1000  seconds  [Davidson  1964].   This  study  aimed 
at  investigating  the  source  and  propagation  mechanisms  by 
means  of  the  variability  of  the  power  spectral  density  as  a 
function  of  1)  the  spectral  distribution  of  the  source  or 
driving  function;  2)  the  coupling  or  impedance  matching  to  the 
medium  for  various  propagating  modes,  and  3)  the  propagation 
characteristics  of  the  medium  including  the  geometrical  effects 
of  the  boundaries . 

The  data  was  taken  in  the  New  Jersey  coastal  plain  area 
(geomagnetic  coordinates  39°  38.5fN  ,  74°  31. 8 'W  and  magnetic 
inclination  71°),  using  a  Rubidium  vapor  magnetometer.   Spectra 
were  then  obtained  by  digital  processing,  using  the  techniques 
of  Blackman  and  Tukey  [1968].   From  these  measurements,  it  was 
found  that  there  are  empirical  expressions  for  the  power  density 
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CnT) 


2 


G(f)  (in  units  on  L-jtz — ]  as  a  function  of  period  as  follows 


for  the  2.5  -  10.0  sec.  period  range  the  expression  is 
G(f)  =  1.26  X  10"1*  T  "   (T  in  seconds). 


For  the  greater  than  10.0  sec.  period  range  the  expression  is 
G(f)  =  6.31  X  10*6  T2-6   (T  in  seconds). 


It  was  also  found  that  these  two  spectral  components  appeared 
to  be  independent  of  time.   In  addition,  these  period  ranges 
contain  spectral  peaks  superimposed  on  the  constant  slopes  at 
5  sec.  and  10  sec.  periods.   A  second  group  of  spectral  peaks 
fall  within  the  17-70  sec.  period  range. 

Later  another  study  made  in  the  same  area  [Herron  196  7] 
reported  that  in  the  1700  -  2300  hour  local  time  frame,  the 
power  spectra  had  the  highest  overall  power  density  levels,  and 
the  strongest  spectral  peaks.   It  was  also  observed  that  the 
log  log  plot  of  Power  vs.  Period  had  a  slope  of  -7.1  dB  per 
octave. 

Recently  the  efforts  of  A.C.  Fraser-Smith  and  J.L.  Buxton 
[1975],  using  a  SQUID  (Superconducting  Quantum  Interference 
Device)  magnetometer.,  enabled  them  to  achieve  higher  sensitivity 
measurements  of  the  geomagnetic  activity  in  the  frequency  range 
.1  -  1H  Hz.   The  measurements  were  taken  at  a  mid-latitude  loca- 
tion, in  Stanford,  California  (geomagnetic  latitude  43.5°N)  on 
the  N-S  field  component  and  spanned  a  2-month  interval  between 
26  January  and  26  March  19  74.   They  found  a  minimum  in  the 
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spectrum  in  the  3  -  7  Hz  frequency  range  while  in  the  .1-3 
Hz  range,  there  is  a  monotonic  decrease  with  frequency.   The 
variation  is  approximately  f~1,25  for  the  .2  -  5  Hz  (pc  1)  range. 
For  the  frequency  range  7  -  14  Hz,  the  magnetic  activity  is 
dominated  by  the  first  Schumann  resonance. 

Finally,  a  report  was  made  on  the  power  spectrum  observa- 
tions made  from  data  collected  in  nine  induction  magnetometer 
observatories  in  widely  separated  locations  about  the  earth 
[Wertz  and  Campbell  1976].   It  was  found  that  the  slope  of  the 
log  (Power  Spectral  Density)  vs.  log  (Period)  was  about  3  at 
low  latitudes  and  4  in  the  Auroral  Zone.   There  is  a  maximum 
of  power  spectral  density  at  auroral  latitudes  and  a  minimum 
at  the  lower  middle  latitudes..   For  10  second  period  measure- 
ments, the  power  spectral  density  varied  from  2  X  10    to 
2  X  10~2  (nT)2/Hz.   For  the  100  second  period  data,  the  power 
spectral  density  varied  from  4  X  10-1  to  2  X  10  2  (nT)2/Hz.   The 
correlation  between  spectral  density  amplitudes  and  sunspot 
number  was  .9  for  12  second  period  data  taken  in  College, 
Alaska.   Earth  surface  conductivity,  ionospheric  variations 
and  magnetospheric  conditions  were  found  to  control  the  ampli- 
tude of  the  signals. 

In  addition  to  these  published  results,  there  is  consider- 
able on-going  work  associated  with  various  Navy  projects,  which 
involves  some  measurements  of  geomagnetic  fluctuations  in  the 
frequency  range  considered  in  this  report.   For  further  details 
see  Ref.  11,  13  and  14. 


19 


III.   EXPERIMENTAL  SETUP 

A.   GENERAL  LAYOUT  OF  EXPERIMENT 

As  mentioned  earlier,  there  are  essentially  two  components 
to  the  geomagnetic  field.   These  are  the  main  field  and  the 
time-varying  components  superimposed  on  the  main  field.   This 
experiment  dealt  with  the  measurement  of  the  time-varying  field 
exclusively,  treating  the  main  field  as  the  stable  or  d.c.  com- 
ponent. 

The  experiment  was  conducted  in  Monterey,  California,  in 
the  community  called  La  Mesa  Village  (Geomagnetic  Latitude 
M-2.5°N).   The  sensors  for  magnetic  field  measurement  were  loca- 
ted in  an  adjacent  wooded  area  subject  to  relatively  few 
magnetic  field  disturbances. 

The  measurement  of  the  field  was  conducted  using  a  Cs  Vapor, 
optically  pumped  magnetometer  with  a  pair  of  sensors  positioned 
a  distance  of  5  5  meters  (180  ft)  apart  on  an  East-West  line. 
The  main  electronics  and  power  supply  equipment  were  housed  in 
Building  318,  61  meters  (200  ft)  from  either  sensor.   Electrical 
power  to  and  signal  from  the  sensors  were  transmitted  via  a 
coaxial  cable  61  meters  long. 

The  geomagnetic  field's  time  variations  were  measured  at 
various  times  both  local  day  and  night,  to  get  a  sampling  of  the 
field  as  a  function  of  time  of  day.   The  general  level  of  geo- 
magnetic activity  ranged  from  very  quiet  to  disturbed  during  the 
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month  of  observations .   Magnetic  activity  index  K£  range  was 
0+  to  8+. 

B.   DESCRIPTION  OF  MAGNETOMETER 

The  sensing  unit  of  the  magnetometer  is  the  Cesium  Vapor 
cell.   The  operation  of  this  sensing  unit  is  based  on  the 
optical  pumping  and  monitoring  technique  [Dehmelt  1957].   This 
technique  has  its  basis  in  the  Zeeman  splitting  of  energy 
levels  of  the  cesium  valence  electron  in  the  presence  of  a 
magnetic  field.   The  energy  levels  are  split  into  sublevels 
which  have  energy  differences  proportional  to  the  total  inten- 
sity of  the  ambient  magnetic  field. 

To  measure  the  energy  differences  between  the  sublevels  we 
note  that  AE  =  hv  and  hence  a  measurement  of  v,  the  frequency 
of  the  emitted  or  absorbed  photons  determines  AE  which  in  turn 
is  proportional  to  the  ambient  magnetic  field  B  since 


AE  =  |g£  B. 


(g  =  Lande  g  factor,  e  =  charge  on  the  electron,  h  =  Planck's 
constant,  m  =  mass  of  the  electron,  B  =  the  ambient  magnetic 
field.) 

The  measurement  of  v  is  accomplished  by  monitoring  the 
optical  transparency  of  the  pumped  Cs  vapor  cell  while  applying 
a  RF  field  slowly  swept  in  frequency  around  v.   When  the  driving 
frequency  corresponds  to  AE,  the  optically  pumped  electrons  are 
redistributed  among  all  of  the  sublevels  leading  to  an  increased 
absorption  of  the  optical  radiation.   Thus  at  this  particular 
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value  of  v,  usually  called  the  Larmor  frequency  vT  ,  we  have 

Li 


AE    =    hvT     =    ff£    B    °^    VT     =77^6 
L        4Trm  L        4irm 


The  constant  of  proportionality  between  the  Larmor  frequ- 
ency and  the  magnetic  field  intensity  for  cesium  is  3.499 
Hz/nT.   The  cesium  sensor  actually  detects  the  magnetic  field 
strength  by  monitoring  the  transmission  of  resonance  radiation 
through  the  optically  pumped  cesium  vapor  cell,  and  produces 
a  signal  whose  frequency  is  proportional  to  the  total  magnetic 
field  intensity  (with  very  fast  response  even  to  large  field 
changes).   Figure  2  illustrates  the  major  components  of  the 
cesium  sensor.   A  more  detailed  description  of  the  optically 
pumped  magnetometer  is  given  in  Bloom  [1961]. 

The  cesium  sensor  employed  throughout  this  experiment  was 
of  the  single  cell  design,  which  has  large  dead  zones.   This 
means  that  the  sensor  can  only  detect  the  magnetic  field  when 
the  sensor  axis  is  within  the  prescribed  cone  in  relation  to 
the  total  field  vector  to  be  measured.   Figure  3  illustrates 
the  active  zones  of  the  single  cell  sensor  and  shows  that  they 
include  the  greater  part  of  one  hemisphere.   Optimum  signal  to 
noise  ratio  is  achieved  when  the  optical  axis  of  the  sensor  is 
45°  from  the  total  field  vector.   Figure  4  shows  the  proper 
alignment  of  a  signal  cell  sensor  in  the  geomagnetic  field, 
for  maximum  signal  to  noise  considerations. 
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Figure  3 
Active  Zones  of  Single  Cell  Sensor 
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Figure  M- 
Sensor  -  Magnetic  Field  Geometry 
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The  sensor  is  effectively  an  oscillator  which  produces  an 
rf  signal  (Larmor  Frequency)  directly  proportional  to  the  total 
field  intensity.   Thus  the  Larmor  Frequency  provides  direct 
information  on  the  field  and  can  be  used  in  conjunction  with  a 
simple  frequency  counter.   However,  the  more  customary  arrange- 
ment is  shown  in  Figure  5 ,  indicating  that  the  Larmor  Frequency 
signal  is  converted  to  a  slowly  varying  voltage,  whose  ampli- 
tude is  proportional  to  the  frequency  excursion  of  the  Larmor 
Frequency,  and  hence  the  magnetic  field  at  the  sensor. 

Finally,  the  following  procedures  had  to  be  followed  for 
reliable  magnetometer  system  operation: 

(1)  The  sensor  was  mounted  in  a  relatively  noise-free  area, 
secluded  from  power  lines,  etc.,  and  positioned  on  a  firm  non- 
magnetic support; 

(2)  The  sensor  was  securely  fastened  to  the  support  to 
preclude  its  movement; 

(3)  The  cable  end  of  the  sensor  was  pointed  in  the  magnetic 
North  direction; 

(4)  For  optimum  signal  to  noise  considerations,  the  sensor 
axis  was  aligned  at  a  45°  angle  from  the  total  field  vector 
(Magnetic  Latitude  42.5°N  and  magnetic  inclination  =  62°  for 
Monterey,  California); 

(5)  The  sensor  was  extended  the  full  six  feet  from  the 
sensor  electronics  box,  and 

(6)  The  sensor  and  sensor  electronics  units  were  sealed  in 
fiberglass  housings  to  prevent  water  intrusion  and  provide  shel- 
ter from  the  wind. 
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Figure  5 
Magnetometer  System 
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The  specifications  supplied  for  the  Varian  Model  49  3  8  mag- 
netometer by  Varian  Associates  are  listed  here: 

Sensitivity  +5  picotesla  (.005  gamma) 

Information  Rate     Continuous  and  essentially  instantaneous 

Range  Continuous  from  20,000-80,000  nT 

C.   DATA  COLLECTION  SYSTEM 

The  major  components  required  to  record  the  data  on  mag- 
netic tape  for  analysis  and  display  at  a  later  time  are: 

(1)  Sensor; 

(2)  Sensor  Electronics; 

(3)  Readout; 

(4)  Differential  Amplifier; 

(5)  Output  Filter,  and 

(6)  Instrumentation  Tape  Recorder. 

Figure  6  shows  the  system  components  required  to  collect 
magnetic  field  data  as  an  analog  signal  (i.e.,  voltage  ampli- 
tude as  a  function  of  field).   In  this  configuration,  it  is 
possible  to  tape  record  the  analog  signal  for  further  analysis. 
Real  time  processing  can  be  carried  on  simultaneously  with  tape 
recording.   A  comparison  of  the  two  spectra  allows  determination 
of  the  tape  recorder  noise. 

In  this  data  acquisition  configuration  called  the  magneto- 
meter mode,  the  sensor  generates  the  Larmor  Frequency  which  is 
proportional  to  the  ambient  magnetic  field.   In  the  Readout, 
the  Larmor  Frequency  is  mixed  with  a  reference  oscillator  in 
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Figure  6 
Data  Collection  System 
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the  mixer.   The  frequency  of  the  reference  oscillator  cor- 
responds to  the  main  field  at  the  given  location.   In  Monterey 
a  reference  oscillator  frequency  of  175,7  80  Hz  corresponding  to 
a  50.3  yT  setting  on  the  Readout  was  used.   The  mixer  signal  is 
a  square  wave  having  a  frequency  equal  to  the  difference 
between  the  Larmor  Frequency  and  the  reference  oscillator 
frequency.   This  square  wave  is  the  input  to  the  discriminator 
driver  and  discriminator,  and  thus  the  change  in  magnetic  field 
is  converted  to  a  d.c.  voltage  whose  amplitude  is  proportional 
to  the  magnetic  field  variation.   In  order  to  use  the  discrim- 
inator output,  a  differential  amplifier  is  used  to  accept  the 
discriminator's  balanced  output  (balanced  output  has  neither 
side  grounded) .   The  output  of  the  differential  amplifier  is 
then  filtered  (Bandlimited  to  .018-300  Hz)  providing  a  suit- 
able signal  for  analog  magnetic  tape  recording.   Thus  this 
system  measures  the  variation  in  the  total  magnetic  field  inten- 
sity. 

The  signal  was  recorded  on  a  Hewlett  Packard  3960  instru- 
mentation tape  recorder  at  a  speed  of  3  3/4  i.p.s.   The  flutter 
compensation  option  was  employed  to  reduce  the  noise  due  to 
vibrations  in  and  near  the  recorder.   The  signal  sent  to  the 
recorder  was  a  bandlimited  (.018-300  Hz)  analog  output  from  the 
discriminator  via  the  differential  amplifier  and  the  filter. 
The  bandpass  filter  had  slopes  on  low  and  high  ends  of  24dB/ 
octave.   All  data  was  recorded  on  Scotch  Brand  instrumentation 
tape,  with  each  side  holding  2  hours  and  2  minutes  of  magnetic 
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field  fluctuation  data.   All  tapes  were  degaussed  prior  to 
using  them  in  data  collection.   The  data  was  recorded  in  the 
FM  recording  mode,  which  improves  amplitude  accuracy,  and  is 
less  sensitive  to  tape  imperfections. 

Some  of  the  pertinent  specifications  for  the  tape  recorder 
are  listed  below: 

Carrier  center  frequency  at  3  3/4  i. p. s.  6.75  kHz 

Pass  band  0  -  1250  Hz 

S/N  ratio  without  flutter  compensation  48 

Distortion  Less  than  1.5% 

d.c.  drift  i.1%  peak  to  peak 

output/ degree 
Celsius 

D.   SCHLUMBERGER  MODEL  1510-03  SPECTRUM  ANALYZER 

This  spectrum  analyzer  is  a  fully-digitalized  high  perform- 
ance instrument  that  provides  real  time  spectrum  component 
analysis  of  analog  or  digital  input  signals  in  256  spectral 
bins  across  one  of  10  selectable  frequency  bands.   The  result 
of  each  of  the  2  56  independent  analyses  performed  by  the  instru- 
ment is  displayed  on  an  internal  C.R.O.  as  a  dot  indicating  the 
power  and  frequency  components  of  the  analysis.   The  bandwidth 
of  each  line  of  resolution  is  dependent  on  the  maximum  frequency 
of  the  selected  analysis  range.   The  1510  translates  time  domain 
signals  into  the  frequency  domain  by  Fourier  transform  techni- 
ques . 
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There  are  five  selectable  input  range  attenuation  switches 
varying  from  .1  V    to  10  V    full  scale.   In  addition,  the 
1510  has  power  spectrum  averaging  selection  controls  to  aver- 
age as  many  as  10  24  individual  spectra. 

The  analyzer  has  an  option  of  two  window  functions  for 
application  to  the  data  in  the  time  domain,  as  an  amplitude 
function  multiplied  by  the  data  samples.   In  effect,  the  entire 
analysis  sequence  is  amplitude  modulated  by  the  window  func- 
tion and  the  result  transferred  to  the  Digital  Fourier  Transform 

analysis  memory.   The  window  function  is  unity  for  the  rectangu- 

n— % 

lar,  and  for  the  Hamming  function  %-%  cos  2tt  YnW  for  eacn  n 

of  the  1  to  1024  data  samples.   The  rectangular  window  has  a 
smaller  equivalent  noise  bandwidth  and  was  thus  employed  in  the 
data  analyzing  scheme. 

For  the  data  analyzed  in  this  experiment,  the  1510  spectrum 
analyzer  was  used  in  one  configuration  only  with  the  following 
settings  kept  constant  for  each  sample: 

Input  Attenuation  .1  V     (OdB) 

Frequency  Range  .1  -  2  5.6  Hz 

Window  Function  Rectangular 

Mode  Continuous 

Spectrum  Axis  X-  log   f ;  Y-dB 

Stop  after  N  depressed 

Averaging  Interval  512 
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E.   CALIBRATION  AND  SYSTEM  TESTS 

A  number  of  system  tests  were  performed  prior  to  data 
collection.   The  first  of  these  tests  involved  the  frequency 
response  of  the  magnetometer  system.   This  test  was  designed  to 
check  the  output  of  the  discriminator  of  the  readout  by  apply- 
ing a  known  input  to  the  sensor  head.   To  accomplish  this  task 
the  sensor  and  readout  were  set  up  as  illustrated  in  Figure  6 . 
Measured  field  changes  were  produced  and  applied  to  the  sensor 
by  a  sinusoidally  stimulated  7  2  turn  coil.   A  sinusoidally 
varying  magnetic  field  component  was  generated  by  an  a.c. 
current  of  known  amplitude  and  frequency  flowing  through  the 
coil  located  as  shown  in  Figure  7. 

The  measured  resistance  of  the  coil  was  1.1  ohms.   The 
inductance  of  the  coil  was  calculated  from  the  equation  for 
inductance  for  a  circular  cross  section  coil. 

L   =  N2ryn[ln(— )-2]  henries 
o       *j    a 

where   N  =  number  of  turns  =  7  2 

r  =  radius  of  coil  =  33  X  10""2m 

a  =  radius  of  circular  cross  section  =  1  X  10  2m 

yQ  =  permeability  of  free  space  =  Utt  X  10  y  g.f 

The  inductance  was  calculated  to  be  7.7  mhenries.   This  set 
up  a  simple  circuit  for  the  stimulation  process  as  shown  in 
Figure  8.   The  impedance  (X  =  R+wL)  is  a  function  of  frequency 
(a>)  and  the  magnetic  field  generated  by  the  coil  is  directly 
proportional  to  the  magnitude  of  the  current  in  the  coil  for 
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Figure    8 

Circuit  Diagram  for 
Magnetic  Field  Generation 
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the  geometry  described  in  Figure  7.   Here  the  generated  field 
lies  along  the  earth's  total  field  vector  so  that  the  coil's 
field  either  adds  or  subtracts  from  the  total  field  vector. 
The  magnetic  field  intensity  along  the  coil  axis  is  given  by: 


NIr2   4tt 
217"  To" 


H  =  g#-  ^  CnT] 


where   N  =  number  of  turns 

I  =  current  in  coil 

r  =  radius  of  coil 

1  =  distance  of  coil  to  sensor 

M-7T 


hit     -  ,  _      .  _  •    amp-turn 

-,  n_2  =  factor  required  for  unit  conversion  from  — *~ 

1U  „,  meter 

to  nT  L 


The  coil  was  "swung"  with  a  sine  wave  of  known  frequency 
and  amplitude  to  change  the  field  a  known  amount.   The  field 
magnitude  change  was  kept  constant  while  varying  the  frequency. 
Figure  9  illustrates  the  results  of  this  test.   Clearly  there 
appears  to  be  a  non-linearity  of  system  response  with  frequency. 
The  test  was  performed  on  two  sensors  with  essentially  identical 
results.   It  was  noted  that  the  frequency  response  was  essen- 
tially flat  from  .1  to  2  Hz,  for  both  sensors  tested.   The 
observed  roll  off  from  2.  to  10.  Hz  was  approximately  1.1  dB 
per  Hertz.   The  exact  reason  for  this  non-linearity  is  not 
known,  but  an  appropriate  correction  for  it  was  applied  to  the 
data. 

The  next  test  involved  determination  of  the  noise  level  of 
the  tape  recorder  that  was  used  in  data  collection.   This  was 
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done  in  two  ways.   First  a  grounded  input  was  applied  to 
channel  No.  1  of  the  tape  recorder  and  the  frequency  spectrum 
was  analyzed  by  the  spectrum  analyzer.   Results  are  found  in 
Figure  10.   The  analysis  was  run  at  tape  recorder  speed  3  3/4 
i.p.s.  through  a  bandpass  filter  set  for  range  .1  to  250  Hz. 
The  spectra  were  averaged  2  56  times.   Major  noise  components 
appeared  at  frequencies  above  the  area  of  interest  and  were 
mostly  power  line  (60  Hz)  components  or  power  line  harmonics. 

Another  test  was  used  to  check  the  reliability  of  the  tape 
recorder  to  reproduce  a  given  analog  magnetic  field  signal. 
This  test  was  designed  to  observe  distortion  of  the  signal  by 
the  tape  recording  process.   Figure  11  gives  the  system  compo- 
nents in  block  diagram  form.   A  real  time  spectrum  analysis 
was  run  while  the  same  signal  was  tape  recorded.   The  tape 
recorded  signal  was  then  played  back  and  also  analyzed  by  the 
spectrum  analyzer.   Graphs  of  the  two  spectra  appear  on  Figure 
12.   The  tape  recorded  spectrum  essentially  duplicated  the  real 
time  spectrum  and  the  tape  recorded  curve  is  within  -2  dB  of 
the  direct  spectrum.   The  small  spike  at  3.3  Hz  was  later  elimi- 
nated by  using  the  flutter  compensation  feature  of  the  tape 
recorder. 

It  is  not  possible  to  put  the  cesium  sensor  into  a  field 
free  region  to  measure  sensor  noise,  since  the  sensor  requires 
at  least  a  20,000  nT  field  to  function.   A  readout  noise  measure- 
ment was  possible  however.   This  measurement  was  made  prior  to 
data  taking,  during  data  taking  and  again  at  its  completion. 
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In  order  to  obtain  equipment  noise  associated  with  the  read- 
out, filter  and  amplifiers,  it  is  necessary  to  provide  a 
constant  frequency  signal  input  as  a  replacement  for  the 
varying  sensor  signal.   This  constant  frequency  signal  must 
be  within  the  range  of  the  total  magnetic  field  present,  so 
that  the  correct  reference  oscillator  is  selected.   With  the 
sensor  disconnected,  the  synthesized  "Larmor  Frequency"  was 
applied  to  the  readout  unit.   The  constant  frequency  synthe- 
sized "Larmor  Signal"  was  beat  against  the  reference  oscillator 
in  the  mixer  and  the  change  in  frequency  is  detected  by  the 
discriminator.   Since  the  change  in  frequency  is  essentially 
zero,  there  should  be  no  signal  emanating  from  the  discrimina- 
tor.  This  test  will  display  any  noise  generated  within  the 
discriminator,  mixer,  filter  and  amplifiers.   Figure  13  shows 
the  results  of  an  averaged  noise  spectrum  taken  over  six  periods 
within  the  span  of  one  month.   The  noise  spectrum  is  thus  quite 
reproducible.   The  low  frequency  end  of  the  noise  spectrum 
exhibited  the  highest  power  content  and  dropped  off  by  20  dB 
at  10  Hz. 

The  power  spectra  obtained  from  the  1510  spectrum  analyzer 
give  only  a  relative  dB  scale  on  the  y  axis.   It  was  necessary 
to  calibrate  the  y  axis  to  the  power  spectral  density  in  y2/Hz 
or  (nT)2/Hz.   This  was  done  following  the  procedure  outlined 
in  the  EMR  1510  Spectrum  Analyzer  Operators  Manual.   The 

log   (nT)2  vs.  log   (frequency)  calibration  is  given  here: 
10^!z~         10 
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Calibration 
Frequency  Range  .1  to  2  5,6  Hz 

Window  Function  (Rectangular)   OdB  correction 
Input  Signal 

(a)  The  calibration  of  the  signal  out  of  the  discriminator 
and  the  differential  amplifier  was  4  mV  per  Hz  of  frequency 
excursion. 

(b)  Filter  gain  was  accomplished  by  an  amplifier  within  the 
filter  and  was  a  factor  of  10.5  X.   Thus  the  signal  was 

10.5/4mV_\   42mV 
\   Hz/     Hz 

(c)  The  cesium  cell  Larmor  Frequency  constant  per  nT  of 
field  intensity  change  is  3.499  Hz  per  nT,  thus 

42mV    3.499Hz    147mV   reference  signal 
Hz        nT  "    nT 

STEP  1;   Relate  input  reference  signal  to  100  mV ,  the 

OdB  attenuation  setting 

147mV   =  InT  +  lOOmV  =  .68:1  =  -3.35dB:   InT 

14  7mV 

Correction   -3 . 35dB 
STEP  2 :   Read  3  correction  factors  from  1510  front  panel 

a.  Frequency  range  .1  to  2  5.6  Hz  -  per  Hertz  cor- 
rection 

Correction   +10dB 

b.  Required  attenuator  setting  to  prevent  overload 

Correction    OdB 
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c.   Window  function  rectangular 

ENBW  =  1  Correction   OdB 

STEP  3:   Add  all  corrections  in  steps  1  to  3  for  total 

correction 

-3.35  dB 
+10.00  dB 

+6.65  dB   TOTAL  CORRECTION 

This  means  that  a  signal  corresponding  to  a  InT   change  in 

field  strength  at  the  1510  input  will  require  correction  of 

6.65dB   at  1510  output.   Thus,  on  the  C.R.O.    "    is  defined 

nz 

to  be  at  the  -6.65dB   calibration  mark.   By  relating  the  input 
signal  to  the  OdB   level  it  is  possible  to  calculate  where  InT 
lies  on  the  y  axis  of  the  dB  calibrated  C.R.O.   The  corrections 
are  necessary  to  take  into  account  what  scales,  window  and 
frequency  range  are  being  used. 
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IV.   EXPERIMENTAL  RESULTS 

A.   INTRODUCTION 

In  a  series  of  preliminary  tests  and  calibration  runs,  as 
outlined  in  the  previous  section,  the  suitability  of  the  mag- 
netometer and  recording  equipment  to  analyze  geomagnetic 
fluctuations  had  been  established.   During  the  month  of  April, 
19  78,  a  total  of  5  3  hours  of  data  was  recorded,  covering 
various  times  of  day  and  a  range  of  geomagnetic  activity  from 
very  quiet  to  moderately  disturbed.   Each  data  run  was  8  5 
minutes  long  from  which  512  individual  spectra  were  computed 
and  then  averaged.   The  width  of  each  frequency  "bin"  was  0.1 
Hz  or  256  bins  between  .1  and  25.6  Hz.   The  averaged  estimate 
of  the  spectrum  of  a  random  input  signal  has  a  Chi-squared 
distribution,  having  two  times  the  selected  number  of  averaged 
spectra  degrees  of  freedom.   The  statistical  characteristics 
result  in  a  95  per  cent  confidence  that  the  estimate  will  be 
within  a  ± . 4dB  accuracy  for  a  512  spectra  average. 

The  observations  made  of  the  magnetic  field  in  Monterey 
were  most  concerned  with  the  0.1  to  10.0  Hz  frequency  range. 
Each  of  the  spectra  taken  the  spectral  power  density  for  this 
range  was  computed  and  displayed.   Comparisons  were  made  of 
the  spectral  power  density  at  one  specific  time  frame  for  a 
number  of  days.   Comparisons  were  also  made  of  spectra  taken 
at  different  times  of  the  day,  and  of  the  spectral  power  den- 
sities as  affected  by  varying  degress  of  field  activity. 

i+6 


B.   TYPICAL  SPECTRA 

Figures  14  through  18  each  represent  the  average  of  512 
spectra  taken  from  a  recording  of  the  earth's  magnetic  field 
at  various  time  periods  during  the  day.   Each  curve  is  a 
typical  spectrum  for  daily  field  variations  during  the  time 
frames  indicated. 

Figures  19  through  21  illustrate  the  spectra  of  five  8  5 
minute  runs  for  local  night  (0000-0200),  local  morning  (0800- 
1000)  and  local  afternoon  (1600-1800),  showing  the  day  to  day 
variations  in  the  spectra.   All  times  are  given  as  Pacific 
Standard  Time. 

Figures  22  through  24  are  the  averages  of  the  five  curves 
given  in  Figures  19,  20  and  21,  except  for  local  night,  where 
only  three  of  the  five  spectra  were  used.   It  was  found  that 
because  one  of  the  curves  had  significantly  more  power  than 
the  other  four,  and  another  of  the  curves  had  significantly 
less  power,  that  it  was  decided  not  to  include  these  two  when 
computing  the  mean.   Otherwise  the  high  power  spectrum  would 
effectively  be  reproduced  as  the  average  spectrum. 

Figure  25  illustrates  the  power  spectra  of  the  1730  to 
19  30  period  for  an  average  day  and  for  a  magnetically  disturbed 
day.   The  disturbed  field  has  lOdB  greater  power  level  in  the 
.1  to  . 6  Hz  range.   The  spectra  nearly  duplicate  each  other 
from  3  Hz  to  10  Hz.   The  Fredricksburg  a  index  had  a  value  of  7 
at  the  disturbed  recording  time  and  a  value  of  3  at  the  time 
the  lower  power  curve  was  taken.   These  index  readings  were 
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provided  by  the  National  Geophysical  and  Solar-Terestial  Data 
Center  in  Boulder,  Colorado. 

C.   COMPOSITE  SPECTRA 

Figure  2  6  is  a  composite  of  the  important  curves  taken 
during  the  course  of  this  experiment.   It  includes  a  five-day 
average  for  local  morning  and  afternoon,  and  a  three-day 
average  for  local  night.   In  addition,  the  equipment  noise 
spectrum  is  displayed.   It  is  evident  from  this  composite  that 
the  local  night  (0000-0200)  spectrum  approached  the  system 
noise  spectrum  most  nearly  of  the  three  periods  observed.   The 
local  night  spectrum  had  significantly  less  power  than  the 
other  two  periods  only  in  the  . 1  to  .  4  Hz  range.   From  .5  Hz 
to  10.0  Hz  all  three  average  spectra  exceeded  the  system  noise 
level  by  at  least  lOdB. 

Local  morning  and  afternoon  had  very  similar  spectra.   At 
most  frequencies  the  two  spectra  were  separated  by  only  3dB 
in  power.   At  all  frequencies  the  morning  spectrum  had  slightly 
more  power  than  did  the  afternoon  spectrum.   Both  spectra  dis- 
played a  definite  -20dB/decade  slope  in  the  0.1  to  1.0  Hz  range, 
Both  curves  also  exhibited  a  drastic  slope  change  from  1.0  Hz 
to  10  Hz.   Here  the  slope  changed  to  positive  at  M-  Hz  and  then 
back  to  negative  at  the  10  Hz  point. 

The  local  night  averaged  curve  had  very  low  power  at  the 
low  frequency  end  of  the  spectrum  compared  to  the  other  two 
periods .   This  was  true  for  the  . 1  to  . 4  Hz  range .   From  . 5  Hz 
to  10.0  Hz,  the  local  night  spectrum  displayed  the  general 
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trend  exhibited  by  the  local  morning  curve.   In  general,  the 
local  night  curve  had  the  lowest  power  density  for  all  frequ- 
encies . 

The  lowest  signal  to  noise  ratio  occurred  at  the  .3  Hz 
point  where  it  had  a  value  of  7.5dB.   A  signal  to  noise  ratio 
of  +10dB  existed  for  all  frequencies  above  .5  Hz. 
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V,   CONCLUSIONS 

A.   DATA  SUMMARY  AND  RECOMMENDATIONS 

The  characteristics  of  the  spectra  agree  in  general  shape 
with  those  illustrated  in  the  Trapped  Radiation  Handbook 
[Cladis,  Davidson  and  Newkirk  1971],  for  the  frequency  range 
0.1  to  10.0  Hz.   The  power  levels  at  these  frequencies  agree 
in  general  with  those  observed  earlier  at  Stanford  [Fraser- 
Smith  and  Buxton  19  75]. 

The  average  spectra  illustrate  the  varying  level  of  magnetic 
fluctuations  as  a  function  of  time  of  day.   For  the  time  periods 
analyzed,  it  was  observed  that  local  night  (0000-0200)  had  the 
lowest  power  density  for  each  frequency  measured.   The  local 
morning  (0800-1000)  spectrum  exhibited  the  greatest  power 
density  for  each  frequency  measured.   The  .1  Hz  point  of  the 
local  morning  spectrum  had  the  greatest  power  density 
(4.  X  10" 2 (nT)2/Hz) .   The  2  Hz  point  of  the  local  night  spectrum 
had  the  least  power  density  recorded  («+.  X  10"s  (nT)  2/Hz)  thus 
spanning  four  orders  of  magnitude. 

It  was  found  that  the  three  average  spectra  converged  as 
frequency  increased  toward  10  Hz,  where  only  2dB  in  power  sepa- 
rated the  night  time  and  morning  power  density  levels. 

During  data  collection,  there  were  periods  when  the  geo- 
magnetic activity  level  was  too  low  to  record.   At  these  times 
only  the  equipment  noise  spectrum  was  observable. 
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The  data  taken  during  the  experiment  and  presented  above 
is  deemed  reliable.   However,  in  order  to  extract  meaningful 
geophysical  conclusions,  it  is  necessary  to  collect  and 
analyze  a  far  greater  amount  of  data  for  different  times  of 
the  year,  different  times  of  the  day  and  for  varying  degrees 
of  geomagnetic  activity  levels.   It  is  therefore  recommended 
that  these  measurements  be  continued  and  expanded,  preferably 
with  the  addition  of  digitizing  equipment  for  computerized 
data  handling  and  analysis. 

B.   EQUIPMENT  IMPROVEMENTS  RECOMMENDED 

The  Varian  49  38  magnetometer  system  provided  reliable  mag- 
netic field  information  in  the  0.1  to  10.0  Hz  frequency  range. 
However,  the  frequency  response  of  the  system  was  not  flat  as 
initially  anticipated,  but  rather  rolled  off  by  8.5dB  from 
2.0  to  10.0  Hz.   The  cause  of  this  roll  off  is  suspected  to 
lie  in  the  mixer-discriminator  circuitry  since  both  sensors 
exhibit  essentially  the  same  performance  properties. 

The  plotted  power  spectra  did  not  show  the  same  fine  struc- 
ture as  observed  on  the  C.R.O.   This  is  primarily  due  to  the 
plotter-analyzer  interface.   Some  fine  structure  elements  were 
lost  in  transmitting  the  analyzer  C.R.O.  display  to  the  plotter, 
due  to  a  response  mismatch  between  the  analyzer  and  plotter. 

The  single  cell  sensor  was  found  to  be  reliable  only  when 
in  a  stable  and  stationary  support.   This  type  of  sensor  can 
not  be  allowed  to  have  large  angular  motion,  and  must  have 
stability  in  its  orientation  to  the  total  field  vector  for 

optimum  signal  to  noise  characteristics. 
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The  tape  recording  of  the  magnetometer  data  allowed  for 
permanent  record  keeping  of  the  data  but  contributed  some 
noise  to  the  system.   The  flutter  compensation  circuit  cancel- 
led out  much  of  the  vibration  noise  and  a  good  reproduction  of 
the  magnetometer  discriminator  output  was  thus  achieved  to 
within  t    2dB.   A  digital  data  collection  scheme  could  greatly 
improve  the  absolute  accuracy  of  the  field  measurements  and 
eliminate  the  discriminator  circuitry  and  its  associated  noise 
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